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Abstract 

We study the hadron resonance gas model and describe the equation of state of QCD and the vacuum expectation 
value of the Polyakov loop in the confined phase, in terms of hadronic states with light quarks in the first case, and 
with exactly one heavy quark in the second case. Comparison with lattice simulations is made. 

Keywords: finite temperature, QCD thermodynamics, heavy quarks, chiral quark models, Polyakov loop 



a; 

CD 



1. Introduction 

Two symmetries are relevant in the study of the ther- 
modynamics of QCD [ 1 ] . On the one hand the chiral 
symmetry is broken at low temperatures, while it be- 
comes restored above a certain temperature T x . In the 
limit of massless quarks, the order parameter for this 
transition is the quark condensate (qq). On the other 
hand, the center symmetry Z(N C ) of the gauge group 
SU(A^ C ) is broken at high temperatures and it is restored 
below T c , known as deconfinement temperature. This 
symmetry controls the confinement/deconfinement of 
color charges |2J]. In the limit of infinitely heavy quark 
masses (gluodynamics), the order parameter for this 

transition is the Polyakov loop L T = {tx c Pe'^ A ° dx °}, 
where Aq is the gluon field and P denotes path order- 
ing. It is accepted nowadays that both transition tem- 
peratures are very close to each other, T x ~ T c , at least 
at zero chemical potential ||3|-|6|]. 

While in the deconfined phase of QCD the quarks 
and gluons are liberated to form a plasma, in the con- 
fined/chiral symmetry broken phase the relevant de- 
grees of freedom are bound states of quarks and gluons, 
i.e. hadrons and possibly glueballs. This means that it 
should be expected that physical quantities in this phase 
admit a representation in terms of hadronic states. This 
is the idea of the hadron resonance gas model (HRGM) 



which describes the equation of state of QCD in terms 
of a free gas of hadrons Jl - 14 ] , 
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with g a the degeneracy factor, g a = ±1 for bosons and 
fermions respectively, and M a the hadron mass. It has 
been presented in 1 3. T^] a similar model to describe 
the Polyakov loop in terms of hadronic resonances with 
exactly one heavy quark. In this communication we will 
ellaborate on these models, and perform a comparison 
with recent lattice simulations. 

2. The hadron resonance gas model and the 
Polyakov loop 

An effective approach to the physics of the phase transi- 
tion is provided by chiral quark models coupled to gluon 
fields in the form of a Polyakov loop Ji, 17 -23t]. Most 
of these works remain within a mean field approxima- 
tion and assume a global Polyakov loop. Based on QCD 
arguments we have shown in 1 15] that a hadronic repre- 
sentation of the Polyakov loop is given by 
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where gi, a are the degeneracies and Aj m = M\ m - nth are 
the masses of hadrons with exactly one heavy quark (the 
mass of the heavy quark itself mi, being subtracted). We 
have shown in II 1611 that Eq. (0 is fulfilled in chiral quark 
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models coupled to the Polyakov, when one goes beyond 
mean field and advocate the local and quantum nature 
of the Polyakov loop. The need of these corrections was 
already stressed in flU [H I2H I2H1 . 

A natural step is to check to what extent the hadronic 
sum rule is fulfilled by experimental states compiled 
in the PDG 0]. The relation © follows in the 
limit nih — > 00, but the heavy quarks in nature have 
a finite mass. Hadrons with a bottom quark would 
be optimal, due to the large quark mass compared to 
Aqcd, but the available data are scarce, so we turn to 
charmed hadrons. Specifically, we consider the lowest- 
lying single-charmed mesons and baryons with u, d, and 
s as the dynamical flavors, with quarks in relative s- 
wave inside the hadron. For mesons, these are usually 
identified with the states (spin-isospin multiplets) D, D s , 
D*(2010) and D*, and for baryons, with A c , 2 C (2455), 
H c , E' c , Q. c , I c (2520), E c (2645), and Q(2770). A total 
of 12 meson states and 42 baryon states. 

The plot in Fig. Q] shows that the lowest-lying states 
fall short to saturate the sum rule, regardless of the 
choice of mass of the charmed quark, m c . This is not 
surprising as any model predicts many excited states on 
top of the lowest-lying ones, as is also the case for light- 
quark hadrons. Adding more states from the PDG does 
not seem practical due to the fragmentary information 
available. Instead we turn in what follows to hadronic 
models. The aim is not so much to have a detailed de- 
scription of the various states but to give a sufficiently 
good overall description of the whole spectrum. 



3. The relativized quark model: trace anomaly and 
Polyakov loop 

Several models in the past have been proposed to de- 
scribe the hadron spectrum and give a prediction for ex- 
cited states not yet included in the PDG. One of them is 
the relativized quark model (RQM) J28, 29]; it is a soft 
QCD model based on a one-gluon exchange at short dis- 
tances and a phenomenological implementation of con- 
finement by a flavor-independent Lorentz-scalar inter- 
action. We will use the spectrum predicted with this 
model to saturate the sum rules. 



3.1. Thermodynamics of QCD 

From the standard thermodynamic relations, the free 
energy, the pressure, and the energy density are given 
by ' 

E T 2 d\ogZ 

F = -pV = -TlogZ, e=- = ^—,(3) 
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Figure 1: Polyakov loop as a function of temperature (in MeV). Lat- 
tice data from (3 for the HISQ/tree action with N^.xN T = 32 3 X 12, 
and from |27] for the continuum extrapolated stout result. The solid 
(brown) line shows the sum rule {2) saturated with the lowest-lying 
charmed hadrons from PDG [ 26]. Dashed lines correspond to the re- 
sult using the RQM spectrum with one b quark and a cut-off A < 
1700 MeV (red line), and A < 5500 MeV (green line). The result 
from the MIT bag model with cut-off A < 5500 MeV is shown as 
a solid (blue) line 1 15]. The results from PDG and RQM have been 
multiplied by a factor L(T) -> e c/T L(T), with C = 25 MeV, which 
comes from an arbitrariness in the renormalization. 



as well as the relation for the trace anomaly 
e - 3p d I p 
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The Hagedorn formula for the trace anomaly follows 
from Eq. ([TJ and the relations given above. It writes 
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where Ki(z) refers to the first order Bessel function. 
n m and rib are the cumulative numbers of mesons and 
baryons (including antibaryons), defined as 



n(M)=Y J g a ®(M-M a ). 



(6) 



is the step function. n(M) represents the number of 
hadrons with mass less than M. Hagedorn proposed that 
the cumulative number of hadrons in QCD is approxi- 
mately given by 



n(M) = A e 



M/T H 



(7) 



where A is a constant and Th is the so called Hagedorn 
temperature. We show in Fig. [2] the cumulative number 
of hadrons with u, d and s quarks computed in the RQM. 
The curves increase up to a cutoff M « 2300 MeV, 
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Figure 2: Cumulative number n as a function of the hadron mass M (in 
MeV), for meson spectrum (blue line), baryon spectrum (red line) and 
total number of hadrons (black line) with u, d and * quarks, computed 
in the relativized quark model of Refs. 12a . 12911 . We also plot a fit of 
the total cumulative number using Eq. (7) (magenta line). 



above which they become approximately flat. The to- 
tal cumulative number can be approximated to the form 
of eq. Q, with A = 0.80, T H = 260 MeV and * 2 /dof = 
0.031, in the regime 500 MeV < M < 2300 MeV. 

As a cross-check of the RQM we can use the spec- 
trum obtained with this model to compute the trace 
anomaly using the HRGM given by Eq. (01. The result 
and its comparison with lattice data is shown in Fig. [3] 
The RQM gives a good description of the trace anomaly 
for T< 180 MeV. 

3.2. Polyakov loop 

The next step is to use the spectrum of hadrons with 
one heavy quark at rest predicted by the RQM to com- 
pute the Polyakov loop using the HRGM of Eq. (0. The 

I with 



total number of hadron states computed in (2S 
one c quark is 1 17 for mesons and 660 for baryons, cor- 
responding to a maximum value of A = M - m c about 
1700 MeV. For hadrons with one b quark, 87 mesonic 
and 776 baryonic states, with a similar upper bound 
for A. In these papers there are some missing states cor- 
responding to baryons of the type csu and bsu, in partic- 
ular H c j,, S' , and Q. c j, baryons. In order to give a pre- 
diction for these states we have used the equal spacing 
rule 13111 . which is based on the approximate relations 



M~. c - M Ai . ~ m s , 
Mq c - My c ~ 2m s . 



M s . -M z 



(8) 



'in order to get a good agreement with lattice data from 13011 . 
a temperature shift of Tq ss 15 MeV is required; ^hrgm(7") = 
^ut6oc(T + Tq). The need of this shift in lattice data was also ob- 
served in [ 30], and it has been attributed to systematic errors in lattice 
coming from extrapolations to the physical light quark masses and the 
continuum limit. 
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Figure 3: Trace anomaly (e - 3p)/T 4 as a function of temperature (in 
MeV). Blue points are lattice data with the asqtad action and N T = 6 
from 1 30]. Dashed lines correspond to a parameterization given in that 
reference of the lattice data with p4 and asqtad actions for N T = 8. 
Brown points are lattice data with the stout action N T = 8 from J2vh - 
The result from the HRGM Eq. {3}, computed with the RQM spectrum 
with u, d and s quarks from Refs. I28tl29ll . is shown as a solid (blue) 
line. A temperature shift of To = 15 MeV to lower temperatures is 
introduced in the lattice data from 133] . 



They follow from the replacement of a light quark (u 
or d quark) in A c by a s quark to get H c , and assum- 
ing m j » m„ , ma- The same replacement in £ c with 
one s quark allows to get H',, and with two s quarks to 
get Q. c . These relations are valid also for baryons with 
one b quark, and they preserve the degeneracy of states. 
After including these missing states, the total number 
of baryons states we consider with one c quark is 1470, 
and with one b quark is 1740. We use an s quark mass 
of 109 MeV (extracted from the lowest-lying hadrons). 

The prediction based on these hadronic states is dis- 
played in Fig. Q] The RQM result is closer to the lat- 
tice data than the naive estimation from the lowest-lying 
hadrons in PDG, but still tend to stay below it, a conse- 
quence of the truncation of states to A < 1700 MeV. In 
order to remove any ambiguities coming from the renor- 
malization prescription used for the Polyakov loop, we 
plot in Fig. 3] the derivative of T log(L(T)) with respect 
to 7". It is noteworthy that the bottom sum rule gives a 
better value than the charm one, as it would be expected 
due to the larger mass of the b quark. 

In order to avoid the cut-off problem in the spectrum, 
we have assumed a raising of cumulative numbers of 
the form n(A) = A A° in the regime 1700 MeV < A < 
5500 MeV, with values A = 5.92 ■ 10~ 9 , a = 3.16 for b- 
mesons, and A = 1.76 • 10~ 24 , a — 8.23 for b-baryons.0 



2 The yellow strip in Fig.f4]is the uncertainty from a combined anal- 
ysis of lattice data from continuum extrapolated stout [ 27], HISQ/tree 
action N T = 12 scale sets r\ and fk, and asqtad scale /j. I14ll . 

3 These numbers follow from a fit of n(A) below 1700 MeV. 
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Figure 4: Comparison of jf(T log(L(T))) (yellow strip) with 
gf(r log( j Ho- Sha? _A,I "' r )) from hadronic states (mesons plus 
baryons): lowest-lying hadrons from PDG (solid brown line, label a), 
RQM states from [28, 29] with quark c (solid red line, label b), and 
with quark b (dashed red line, label c), and MIT bag model estimate 
including states up to A = 5500 MeV (solid blue line, label d) fiUl . 



The analytic formula to be applied in this regime is 



L(T) 



5»(A) e _ A/r 
<9A 



(9) 



which gives a contribution to be added to the numer- 
ical one below 1700 MeV. The result is displayed in 
Fig. Q] with a dashed (green) line. The estimated ef- 
fect of adding the states up to A = 5500 MeV in the 
RQM leads to a result quite consistent with the one ob- 
tained with the MIT bag model discussed in 11511 . with 
the same cut-off. We have also checked that this result 
is numerically indistinguishable from that obtained by 
extending the analytic raising up to A — > oo (no cut-off). 

4. Conclusions 

The hadron resonance gas model was proposed as a 
simple model to describe the confined phase of QCD 
in terms of a free gas of hadronic states (mesons and 
baryons). Using the hadronic spectrum obtained by the 
relativized quark model of Ref. li28L l29h . we get a good 
description of lattice data for the trace anomaly up to 
T = 180 MeV. A different version of the hadron reso- 
nance gas model has been proposed in 115 to de- 
scribe the vacuum expectation value of the Polyakov 
loop in terms of hadronic states with exactly one heavy 
quark at rest and several dynamical quarks. The lowest- 
lying charmed mesons and baryons included in the PDG 
gives for the Polyakov loop a value well below lattice 
data, and this suggest the need for inclusion of more 
states. When using the spectrum predicted by the rel- 
ativized quark model and the MIT bag models up to a 



cut-off A = 5500 MeV, we get a good description of 
lattice data for the Polyakov loop in the confined phase. 

This work opens the possibility of a Polyakov loop 
spectroscopy, i.e. using the Polyakov loop in funda- 
mental and higher representations to deduce multiquark 
states, gluelumps, etc, containing one or several heavy 
quark states. 
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